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ABSTRACT 

We have succeeded in obtaining magnetized equilibrium states with differential rotation and 
differential toroidal magnetic fields. If an internal toroidal field of a protoneutron star is wound 
up from the initial poloidal magnetic field by differential rotation, the distribution of the 
toroidal magnetic field is determined by the profile of this differential rotation. However, 
the distributions of the toroidal fields in all previous magnetized equilibrium studies do not 
represent the magnetic winding by the differential rotation of the star. In this paper, we inves¬ 
tigate a formulation of a differential toroidal magnetic field that represents the magnetic field 
wound up by differential rotation. We have developed two functional forms of differential 
toroidal fields which correspond to a u-constant and a j-constant field in analogy to differ¬ 
ential rotations. As the degree of the differential becomes very high, the toroidal magnetic 
field becomes highly localized and concentrated near the rotational axis. Such a differential 
toroidal magnetic field would suppress the low-T/|lV| instability more efficiently even if the 
total magnetic field energy is much smaller than that of a non-differential toroidal magnetic 
field. 
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1 INTRODUCTION 

A rapidly rotating protoneutron star are produced by a core collapse 
supernova, accretion-induced collapse white dwarf. T hey can be 
subje ct to well-known non-axisymmetric instabilities jAnderssonl 
I2OO3I) . Stellar deformations caused by rotation and the non- 
axisymmetric instabilities can generate strong gravitational wave 
signals, with the global m — 2 instabilities particularly related to 
gravitational wave production. Thus, a rapidly rotating protoneu¬ 
tron star is a good source of gravitational radiation. 

The development of non-axisymmetric instabilities is charac¬ 
terized by the energy ratio /3 = r/| VK|, where T is the rotational 
energy and W is the gravitational potential energy of the rotating 
star. Studies of these inst abilities were obtained in the nineteenth 
century (e.g.. lBrv^l 18891) . Uniformly rotating incompressible flu¬ 
ids exhibit the m = 2 f-mode dynamical instability when the value 
of 0 satisfies the threshold j3 > 0.27 (see also IChandrasekhaj 
1 19691) . On the other hand, there exists a secular instability which 
is driven by dissipative processes, such as viscosity and gravita¬ 
tional radiation, which grows when the value of P is larger tha n 
0.14 (see also numerical calculations bv IShibata & Karinoll2004l) . 
This threshold is easier to satisfy than the threshold for dynami¬ 
cal bar-mode instability, but the growth time-scale is much longer 
because the typical time-scale of the diss ipation process is m uch 
larger than that of a dynamical time-scale dUai & Shanirdl 19951) . 
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The instability situation changes if we introduce sufficiently 
large differ ential rotation. The linear analysis of a rapidly rotat¬ 
ing star by iKarino & Eriguchiil ( l2003h showed that the threshold 
of the dynamical instability decreases when the degree of the dif¬ 
ferential rotation is sufficiently large. Many numerical simulations 
have confirmed that m = 1 or m = 2 non-axisymmetric insta¬ 
bilities can appear at a much lower {3 than 0.27 if a sufficiently 
large differential rot ation is present (for example, fi ~ 0.14 in 
ICentrella et al.l200i]) . This non-axisymmetric instability is referred 
to as low-r/|W| instabili t y, and it can occu r even for /3 ~ 0.01 
dShibata et alj|200^ l2003l : Isaiio et aDl2003l) . Although numerical 
simulations have shown their existence, the physical origin of low- 
T/|IU| instabilities remains unclear. The instabilities are associated 
with the presence of a corotation resona nce point inside the star 
dWatts et al.ll2005l : rSaiio & Yoshidall2006h . The wave pattern speed 
of an unstable mode equals the local background fluid rotational ve¬ 
locity at the point. Corotation resonance has long been known to be 
a key factor for instabilities in various astrophysical contexts, such 
as the corotational instabiliti es for thin accretion discs, the spiral 
pattern in galaxies dShull992h. and the Papaloizou-Pringle instabil - 
ity for accretion tori dPanaloizou & Pringle! 1984 | : IFu & Laj2011ah . 
By performing ti me evolutions of linear perturbati ons of differen¬ 
tial rotating stars, IPassamonti & AnderssonI d2015l) recently found 
that the I = m = 2 f-mode becomes unstable as soon as a corota¬ 
tion point develops inside the star. The low-T /\W\ instabilities are 
thus a subclass of local shear instabilities. 

The effect of magnetic fields is a significant issue concerning 
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2 K. Fujisawa 

these instabilities. Many numerical simulations have shown that 
low-r/|VK| instability may develop in some pro to-neutron stars 
with rapid rotation during stella r core-collapses iOtt et all 1200.51 
bOOTtlShibata & Sekiguchii2005h because of their sufficiently large 
differential rotations. On the other hand, recent magnetohydrody- 
namical simulations have shown that the dyna mical bar-mode insta¬ 
bility can be supp ressed by magnetic force jCamarda et al.|[200^ : 
iFranci et al.ll2013h . Since large toroidal fields can be wound up 
from an initial poloidal magnetic field by di fferential rotation and 


Shibata et alJl2006l: Obersau 

inser et al.ll2006l:lobercaulinBer et al.l 

20091: Masada etalJ 20121: 

Sawai et al.l 2013b. ranidly rotating 


proto-neutron stars are expected to have large toroidal magnetic 
fields, which would suppress the low-T /\W\ instability. 

Contrasting with the many numerical simulation studies, 
IFu & Lail jlOllbh investigated the effect of a purely toroidal mag¬ 
netic field on the low-T /\W\ instability by applying linear analysis 
to a c ylindrical star with a simple model of the toroidal magnetic 
field JOstrikeill 19651 : [^aiio & Yoshidall2006h . They found that the 
toroidal magnetic field suppresses the low-T /\W\ instability when 
the total toroidal magnetic energy M is in the order of 0.2 T or 
larger. Such magnetic energy corresponds to toroidal fields of a few 
10^® G in their magnetic field model. This study also found that the 
corotation point is split into two slow magnetosonic resonances in 
the presence of strong toroidal fields. Therefore, the unstable mode 
in the presence of the toroidal field is determined by both corotation 
and slow magnetosonic resonances inside the magnetized star. 

Very recently, however, iMuhlberger et al.l ^2014) performed 
accurate general relativistic magnetohydrodynamical simulations 
and found that the low-T/|VF| instability woul d be suppressed b y 
a weaker magnetic field energy than given in IFu &Lail l l2011bh . 
In their simulation models, the large toroidal magnetic field (~ 
10^®G) is wound up from the initial poloidal field (~ lO^'^G) by 
differential rotation. Since the toroidal magnetic fields are wound 
up in their models, the distributions of the toroidal magnetic field 
are determined by the profile of the differential rotation. The distri¬ 
butions of the toroidal magnetic fields are different from the simple 
model used in iFu & ( 1201 Ibh . The maximum strength of the lo¬ 

cal toroidal magnetic field reached approximately 10^^ G in their 
simulations, but the total magnetic energy peaked at 0.56% of the 
kinetic energy T. Although thi s energy ratio (A4 /T ~ 5.6 x 10“®) 
is much smaller than that in IFu & i3 ( l2011bl) (MjT ~ 0.2), 
the low-T/1FFI instability was stabilized by the toroidal magnetic 
fields in their simulation models. They also showed that slow mag¬ 
netosonic resonances appear in their simulation models (fig. 10 in 
their paper). This result shows that both the total magnetic energy 
and the structure of the toroidal field are important for suppres¬ 
sion by magnetic fields because low-T/|IF| instabilities signifi¬ 
cantly depend on the value s of the local rotational velocity (degre e 
of the differential rotation: Ishibata et alJl2003l) . lFu & Lm | |201 Ibh . 
however, calculated only one simple toroidal magnetic field model 
which did not represent the wound magnetic field. Therefore, we 
must consider various toroidal magnetic fields in order to study the 
linear analysis and address the suppression systematically. 

Although many equilibrium states with purely toroidal mag- 


(e.e..lMiketinadll973l:lLander & Jonesll200S 

:|Fu & Lail201 id) and 

general relativistic (e.g.. 

Kiuchi & Yoshida 

I 2 OO 8 I: Yasutake et al.l 

120101: iFrieben & Rezzolla 

2 OI 2 I: Pili et alj 2014) frameworks, all 


of their toroidal magnetic field configurations are essentially same- 
class solutions. When we calculate a stationary axisymmetric and 
barotropic rotating star with purely toroidal magnetic fields, we 


obtain two arbitrary functions from the integrability condition of 
the magnetohydrodynamic equation that are related to the rota¬ 
tion law of the angular velocity and the distribution of the toroidal 
magnetic field. Rotating equilibrium states have been studied and 


many rotation laws have been investi gated in both Newtonian 
and general rela ti vistic frameworks (e.g. jEriguchi & MuelleJ 19851 : 
iHachisij Il986al : I Komatsu et alj Il989d : Galeazzi et al. 2012h . In 
contrast, only a few works have studied purely toroidal magne¬ 
tized equilibrium states. Almo st all previous works have employed 
only a simple functional for m jMiketinaj[l973|; |lGuchi_&Yqshida| 


2008l:lLander & Jonesl200^:lYasutake et alj201(]l:rFu & Lail201 id: 


Frieben & Rezzollal2012l : Pili et alj2014l) . The functional form of a 
toroidal magnetic field is considered rigid in analogy to the rotation 
law. Although a rigid toroidal field is simple and widely accepted, 
the distributions of the toroidal field do not represent those wound 
up by differential rotation. Thus, we need to construct a differen¬ 
tial toroidal magnetic field which represents magnetic winding by 
differential rotation and to treat the linear analysis and the suppres¬ 
sion more systematically. As a first step, this paper investigate new 
functional forms of a differential field and obtains new solutions 
of purely toroidal magnetized equilibria. These solutions will be 
useful for equilibrium models for linear analysis. 

The remainder of this paper is organized as follows. In Section 
2, we describe the new functional forms of the arbitrary functions. 
In Section 3, we present results from numerical calculations. Dis¬ 
cussion and conclusions are given in Section 4. 


2 BASIC EQUATIONS AND FORMULATIONS 

We em plo yed the formulation d eveloped by iKiuchi & Yoshidal 
( l2008h and luander & Jon^ l l2009ll . We will first briefly summarize 
the formulation below while the remainder of the section focuses 
on explaining the new functional forms in detail. Both cylindrical 
coordinates {vj, tp, z) and spherical coordinates (r, 9, ip) are uti¬ 
lized in this paper. 


2.1 Assumptions and basic equations 

As per previous works, we make the following assumptions for 
magnetized stars. 

(i) The system is in a stationary state, i.e. d/dt = 0. 

(ii) The rotational axis and magnetic axis coincide. 

(iii) The configuration is axisymmetric about the rotational axis, 
i.e., d/d<p = 0. 

(iv) There are no meridional flows or poloidal magnetic fields, 
i.e. the magnetic field is purely toroidal. 

(v) The star is self-gravitating. 

(vi) The system is treated in the framework of non-relativistic 
physics. 

(vii) No electric current is assumed in the vacuum region. 

(viii) There is no current sheet, i.e. the magnetic field continues 

smoothly. 

(ix) The barotropic equation of state is assumed. 

These assumptions are adopted here for simplicity. For a barotropic 
equation of state, we use a polytropic equation of state 

P = (1) 

where P and p are the gas pressure and the density, respectively, 
and K and N are constants. 
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Magnetized stars with a differential field 3 


Under these assumptions, the equations that describe this sys¬ 
tem are the Euler equations 

—SJp = —^ 4^g -\- F (2) 

P 

F — — + 0 -\- 0 Gz, (3) 

'UJ 

the integral form of the Poisson equation 

Mr) = -G f (5) 

J \r-r'\ 

and stationary Maxwell’s equations, 

V • B = 0, ( 6 ) 


V X B = drr —, 
c 


(7) 


V X B = 0, (8) 

where (j>g, Vg,, j, c, G and E are, respectively, the gravita¬ 
tional potential, the toroidal velocity, the current density, the speed 
of light, the toroidal magnetic field, the gravitational constant, and 
the electric field. Using these e quations, we obtain th e first integral 
of the Euler equation (see e.g. jLander & Jon^l2009h as 


/ 


dp 

P 


— 4^9 “h 4r — 4m -I“ c, 


(9) 


4r 



LO duJ , 


( 10 ) 


2.2 Functional forms of arbitrary functions and their 
physical meanings 

First, we will review functional forms of Q and discuss the rotation 
laws. Th e following three types of rotation law ar e usually used 
(see e.g.. lEriguchi & MuelleJl98^ : lHachisulll986ah : 


(i) rigid rotation: 

n'^{zu) = ffo. (14) 


(ii) w-constant rotation: 


rd{vj) = 


Vo 


+ d'^)’ 

(iii) y-constant rotation: 


(zu^ + d?f ’ 


(15) 


(16) 


where flo, Vq, jo, and d are constants, d represents the degree of the 
differential rotation. As d increases for the latter two laws, their dif¬ 
ferential rotations approach rigid rotation. As d decreases towards 
zero, the latter two rotation laws approach the true rotation denoted 
by their respective names. That is, Vg, becomes strictly constant in 
space for the second law (u-constant rotation) and the specific an¬ 
gular momentum per unit mass, j = zu’^td, becomes exactly con¬ 
stant in space for the third law (j'-constant rotation). For these cases, 
we can obtain the explicit forms of the rotational potentials as fol¬ 
lows. 


(i) For rigid rotation: 

4r{u7) = , (17) 

(ii) For n-constant rotation: 


where C is an integral constant, and (ju and (j)M are the rotational 
potential and the magnetic potential respectively. Q,{vj) and /(A) 
are arbitrary functions of vo and A, respectively, and A is a scalar 
function defined by A = p-uzd. These arbitrary functions are related 
to the toroidal velocity and the toroidal magnetic field according to 

Vg, = zaj^I{-cu), ( 12 ) 




7(A) 

'UJ 


(13) 


Therefore, there appears to be two arbitrary functions when we cal¬ 
culate a rotating star with a purely toroidal magnetic field. Note that 
we can choose any kind of functional forms of these potentials re¬ 
spectively under the assumption of no poloidal magnetic field and 
no meridional flow, because the toroidal magnetic field does not 
produce a force that acts directly onto the rotation and the rota¬ 
tion does not change the toroidal magnetic field. We can obtain one 
rotating magnetized equilibrium solution by fixing these two func¬ 
tional forms. The arbitrary function Q corresponds to the rotation 
law of the rotating star. We can obtain various differential rotat¬ 
ing stars by changing the functional form of fl. On the other hand, 
the physical meaning of the arbitrary function I is less clear. In 
all previous works, a simple functional form was adopted with no 
investigation of new functional forms that express various toroidal 
magnetic fields such as differential rotations. This paper discusses 
new functional forms of I that correspond to the toroidal magnetic 
field wound up by the differential rotation. 


4r{'^) = 


Vq ln(ti7^ -I- d^) 


(iii) For j-constant rotation: 


4r{w) = — 


jo 


2(ci7^ -I- d^) 


(18) 


(19) 


These types of rotation laws have also been studi ed in general 
relativistic framewor ks (e.g.. iKomatsu et H] Il989j f9). Recently, 
lOaleazzi et al.l ( 1201 2 h studied rotation laws systematically and ob¬ 
tained a new functional form. 


fl('ro) = fic 



( 20 ) 


where flc and 7 are constants. We ca n invoke a wide rang e of rota¬ 
tion laws by changing the value of 7 dOaleazzi et al .12013) . Almost 
all works concerning the low-T /\W\ instability have applied the j- 
const ant functional form for di fferential rotations of initial models 
(e.g.. lShibata et alj|200ll2003h . 

In contrast, only the following simple power-law type 
functional form of I has been adopted in p revious works 
dKiuchi & Yoshid^ l2008l: lUander & Jones! l2009l : lYasutake et al.l 
l2010l:lFrieben & Rezzollall2012l:|Pili et alj20l4 


/(A) = loX" 


( 21 ) 


where fc is a constant. Since we have assumed that the magnetic 
field continues smoothly, the regularity of the toroidal magnetic 
field is required. The parameter k must be fc > 1, because the 
functional form with k > 1 satisfies the regularity condition at the 
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centre —>■ 0, to —>■ 0). We can derive the magnetic potential 

using this power-law functional form as 






II k 

47r 2fc - 1 

II k 

47r 2fc - 1 


^2.-1 

/ 2\2k-\ 

[pzu ) 


( 22 ) 


As we have seen, the exponents of the rotational potential are ~ 
TO^ (rigid), <j)R ~ to“^ (v-constant), and (hi} ~ to“^ (j'-constant). 
The exponents of rotational potentials of the differential rotations 
are negative, while that of rigid rotation is positive. The exponent of 
the rotational potential characterize the rotation law. On the other 
hand, the magnetic potent ial becomes (I)m ^ > 1) when 

the star is incompressible jKiuchi & Yoshidal2008h . Since k is lim¬ 
ited to fc > 1 by the boundary condition, the exponent of this mag¬ 
netic potential is always positive. Therefore, we consider this func¬ 
tional form of magnetic potential as a rigid toroidal magnetic field 
in analogy to the rotation law. 

If the toroidal magnetic field is wound up from the initial 
poloidal fields by the differential rotations, the distribution of the 
toroidal magnetic field is determined by the profile of this differ¬ 
ential rotation. If we neglect the meridional flow, the evolution of 
the toroidal magnetic field is described by the following induction 
equation: 

^ = {V X (t; X B)}^ ~ ( 23 ) 

If the initial magnetic field is purely poloidal (Bp) and Bp only has 
a constant TO-component (Bp = Bq e.^ -|-0 -|-0 e^), the poloidal 

field is wound up by the differential rotation and the profiles of the 
toroidal magnetic field become 




Bo —(toD) 
aw 


TO ^ (u—const.), 

to“^ (j— const.). 


(24) 


Therefore, the toroidal magnetic field wound up by the differential 
rotation must satisfy the magnetic winding in equation d24b . but 
we have not obtained such toroidal magnetic field distributions by 
using a rigid toroidal field (equation ED. Since the rotation and 
the toroidal magnetic field are not necessarily consistent with each 
other under the assumption of no meridional flow and no poloidal 
magnetic field, we can select the functional forms of /(A) that are 
most compatible with the chosen rotational potential. 

Here, we introduce two new functional forms of I whose 
toroidal magnetic fields satisfy equation We call the func¬ 
tional forms differential toroidal fields in analogy to rotation law 
and rotational potentials. The functional form of a type 1 differen¬ 
tial toroidal field is 

and the functional form of a type 2 differential toroidal field is 


where is a constant, h represents the degree of the differential 
toroidal magnetic field. The magnetic potentials of these functions 
are 


<)>m(A) 


/q ^ ln(A -t- 

4^ \ 2 


1 

2(A + /i2) /’ 


for a type 1 differential magnetic potential and 


<)>m(A) 


47r 2 (A-b/i2)2 ’ 


(27) 


(28) 


for a type 2 differential magnetic potential. Since the exponents 
of these differential magnetic potentials are ~ to~^ in type 1 and 
~ to”"* in type 2, we can obtain a wide range of the distributions of 
the toroidal magnetic field using these functional forms, including 
distributions which meet the condition of equation l l24b (see the 
profiles of toroidal magnetic fields in fig.[D- 


2.3 Numerical settings and accuracy check 

Following previous works jKiuchi & Yoshid j 1200 Sl : 
iLander & JonesI 20091), we adopted Hachisu’s self-consistent 
scheme jHachisiJ 1986al : lHachisull 1986bh . which is a very powerful 
numerical scheme for obtaining a rotating star in equilibrium. 
In this scheme, we fix the axis ratio q = Tpoi./r^ to obtain 
the equilibrium configuration, where rpoi. and are the polar 
radius and the equatorial radius, respectively. When we calculate 
magnetized rotating equilibria, we need to fix another parameter. 
We choose Iq in the arbitrary functions as the parameter in this 
paper. Thus, we fixed both q and Iq during our numerical iteration 
cycles and obtained one equilibrium state after the iteration. 

The physical quantities are transformed to be dimensionless 
in actua l numerical computat i ons. W e adopted the dimensionless 
forms in iFuiisawa & Eriguchil j2013h by utlising the maximum den¬ 
sity pmax, the maximum pressure Pmax, and the equatorial radius 
r, as follows: 




(29) 


a ‘iiiGpl,., 



Pmax 


(30) 


where a is introduced to ensure a unity distance from the cen¬ 
tre to the equatorial surface of the star. Other physical quantities 
are also transformed to b e dimensionless (see the Appendix in 
iFuiisawa & Eriguchill2013h . A term with a hat 0 denotes a dimen¬ 
sionless quantity. 

To see the global characteristics of magnetized rotating equi¬ 
libria, we define some integrated quantities as follows: 



(31) 



(32) 



(33) 


U = ivn, 


(34) 


M = ^ j BldV, (35) 

where W, T, H, U, and M are the gravitational energy, the rota¬ 
tional energy, the total pressure, the internal energy, and the mag¬ 
netic energy res pectively. We also defi ne the volume averaged mag¬ 
netic field Bave jFuiisawa et alj20T^ as 

_ ijB^pdV 

- Y jdv ~ V 

where V is the volume of the star. This term is useful when we 
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evaluate the localization of the magnetic field. We also calculate 
the toroidal Alfven velocity as 


^4 


Btp 


(37) 


Note that we checked the numerical convergence by using the Virial 
theorem (Appendix A) and ensured a sufficient number of mesh 
points for numerical computations in this paper. 


3 NUMERICAL RESULTS 

This section shows our numerical solutions with differential rota¬ 
tion and a dijferential toroidal field. We fixed the polytropic index 
= 1 as a simple equation of state for a proto-neutron star. We 
adopted the y-constant rotation law with d = 0.3 as the differen¬ 
tial rotation. The low-r/|VU| instability of this differential rotat- 
ing star commences when the value of P becomes larger than 0.03 
dShibata et al.ll2003h . We fixed the axis ratio q = 0.8 and obtained 
the rotating equilibrium state whose /3 is 0.04. Panel (a) in fig. 
is a density contour of this rotating star without toroidal magnetic 
fields. The isopycnic surfaces are distorted by the strong rotation. 
This rotating star gains the low-T /\W\ instability if it does not have 
a toroidal magnetic field. 

To show the profiles of the re-dimensionalized toroidal mag¬ 
netic fields, we used a neutron star’s typical values of mass (M = 
IAMq) and central density (pmax = 10^®g/cm®). The maximum 
strength of the magnetic field is approximately a few times 10^® 
G and the maximum value of the Alfven velocity is a few times 
10®cm/s in this model. 


3.1 Dijferential toroidal magnetic field 

First, we will compare one non-magnetized equilibrium state and 
six magnetized equilibrium states with rigid or dijferential toroidal 
fields. The numerical models and results are tabulated in table [T] 
and plotted in figures [T] and The relative value of the Virial rela¬ 
tion (VC) in all solutions is very small (see accuracy verification in 
Appendix A). We calculated many magnetized equilibrium states 
and displayed some solutions whose ratios of the rotational energy 
to magnetic energy are equal (M/T — 1.0 x 10“^). This energy 
ratio is smaller th an Fu & Lai’s criterion (M jT ~ 2.0 x 10“^), but 
larger than that of lMuhlberger et al] (1201 4h {M /T ~ 5.6 x 10“®). 
These toroidal magnetic fields could suppress the low-T/|FF| in¬ 
stability. 

Since the values of /3 for these solutions are almost the same, 
their magnetic energies are also almost the same. As seen in Table 
[T] however, the ratio Bmax/73ave is different for each model. The 
ratio for the rigid toroidal solution with fc = 1 is 1.88. This value 
is almost as the same as those for the type 1 differential toroidal 
solution with h = 0.1 (1.70) or the type 2 differential toroidal 
solutions with h = 0.1 (1.8) and h = 0.2 (1.99). On the other 
hand, the ratios of the rigid toroidal solution with k = 2 and the 
type 2 differential toroidal solution with h = 0.3 are 2.48 and 2.11, 
respectively. These larger values indicate that the toroidal fields of 
these two models are more localized than those of other solutions. 

The profiles of the toroidal magnetic fields are different 
and rely heavily on their type. Fig. [T] displays the profiles of 
dutid / d-uj (left-hand panel) and the re-dimensionalized B^p (right- 
hand panel). The curves represent the following: (i) the thick solid 
curve shows a type 2 differential field with h — 0.1, (ii) the thick 
dashed curve shows a type 2 differential field with h = 0.2, (iii) the 


thick dotted curve shows a type 1 differential field with h = Q.l, 
(iv) the thin solid curve shows a rigid field with k = 1 and (v) 
the thin dashed curve shows a rigid field with k = 2. The mag¬ 
netic winding by the differential rotation {ddjClIdj in the left-hand 
panel) reaches a maximum value at the rotational axis ut = 0 and 
the value decreases rapidly near the centre. On the other hand, the 
profiles of the rigid toroidal fields reach maximum values in the 
middle regions [ui ~ 0.5 with A: = 0 (curve iv) and ro ~ 0.6 with 
k = 1 (curve v)]. These distributions do not represent the mag¬ 
netic winding by the differential rotation {dvo^jd-uj in the left- 
hand panel). The profile of the type 1 differential toroidal field is 
gradual. The ratio of Bmax/Sav. is the smallest among these seven 
solutions. This functional form of I corresponds to -constant 
toroidal fields such as the u-constant rotation law, because the func¬ 
tional form of the type 1 toroidal field is 7 ~ zu{Btp = const.) 
when the star is incompressible. 

In contrast, the value of the type 2 differential toroidal fields 
reach a maximum near the rotational axis [to ~ 0.2 with h = Q.2 
(curve ii) and to ~ 0.1 with h = 0.1 (curve i)]. As the value of h 
decreases, the degree of the differential toroidal field becomes large 
and the maximum point reaches the rotational axis. These type 2 
differential toroidal fields resemble the magnetic winding by the 
differential rotation. In particular, the toroidal field with h = 0.1 
closely matches the magnetic winding (fig. [TJ. Therefore, we can 
obtain the toroidal field wound up by differential rotation by using 
type 2 differential toroidal fields under the assumption that initial 
profile of the magnetic field is uniform (equation |24j . 

Fig-IHshows the contour maps of p and B^ for these solutions. 
Panel (a) displays the density distribution of the non-magnetized 
solution. Since the magnetic energy of our solutions is smaller than 
the gravitational energy (A1/|FF| — MjT y. T/|1U| ~ 10“"^), 
all solutions in table [T] have similar isopycnic surfaces. The other 
panels show distributions of the toroidal fields. Panel (b) shows 
the contour of a rigid toroidal field with A: = 1. This configu- 
ration is the same as those in previous works (e.g., see fig. 2 in 
lUander & Jonesll2009l) . Panels (c) and (d) display a type 1 differen¬ 
tial field with h = 0.1 and a type 2 differential field with h = 0.1, 
respectively. As we see in fig. [T] the toroidal magnetic fields of a 
type 2 differential field are concentrated and localized near the ro¬ 
tational axis. 


3.2 Highly localized differential toroidal magnetic field 

Next, we calculated the type 2 differential toroidal field by chang¬ 
ing the value of h and compared the resulting models and the rigid 
toroidal field model. We show one rigid toroidal model and six dif¬ 
ferential toroidal models. Since h determines the degree of the dif¬ 
ferential toroidal magnetic field, the toroidal magnetic field is lo¬ 
calized and concentrated near the centre when the value decreases. 
We calculated many equilibrium states by changing the value of h 
and 7o, and obtained one solution sequence. The numerical results 
of the solution sequence are tabulated in Table and the profiles 
of the re-dimensionalized toroidal magnetic fields and the Alfven 
velocity are plotted in fig. In this sequence, the local maximum 
values of the differential toroidal magnetic fields are comparable 
to those of a rigid field at that point (see fig. [3}. 

As seen in fig. the location of the maximum point ap¬ 
proaches the centre as the value of h decreases. This is because 
the degree of the differential toroidal magnetic field increases. The 
maximum point of the rigid toroidal field is located at to ~ 0.5 
and that point in the differential toroidal fields becomes to ~ 0.4 
{h = 0.5), TO ~ 0.2 (h = 0.2), to ~ 0.1 {h = 0.1) and 
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Table 1. Numerical results for each model. All solutions have almost the same magnetic energy. 


q 

^0 

k or h 

P 

MjT 

-Bmax(G) 

-^max/ -Save 

VAmax (cm/s) 

VC 

No magnetic 

0.8 

0 

- 

4.00E-2 

O.OOEO 

0 

- 

- 

l.llE-5 

Rigid toroidal 

0.8 

7.12E-3 

k = l 

4.02E-2 

l.OOE-2 

3.67E16 

1.88 

4.14E8 

1.12E-5 

0.8 

2.95E-1 

k = 2 

4.02E-2 

l.OOE-2 

4.84E16 

2.48 

5.78E8 

1.12E-5 

Type 1 differential toroidal 

0.8 

9.27E-4 

h = 0.1 

4.02E-2 

l.OOE-2 

3.30E16 

1.70 

3.15E8 

1.12E-5 

Type 2 differential toroidal 

0.8 

3.32E-4 

h = 0.3 

4.02E-2 

l.OOE-2 

3.50E16 

1.80 

3.20E8 

1.12E-5 

0.8 

1.72E-4 

h = 0.2 

4.02E-2 

l.OOE-2 

3.85E16 

1.99 

3.45E8 

1.12E-5 

0.8 

8.30E-5 

h = 0.1 

4.03E-2 

l.OOE-2 

5.39E16 

2.77 

4.80E8 

1.15E-5 



0 0.2 0.4 0.6 0.8 1 


A 

G3 



Figure 1. Left: the prohle of dvo^ldvo with d = 0.3 j-constant rotation. Right: the profiles of on the equatorial plain. The vertical axis denotes the 
re-dimensionalized toroidal magnetic field in units of Gauss, (i) Type 2 differential toroidal field with h = 0.1. (ii) Type 2 differential toroidal field with 
h = 0.2. (iii) Type 1 differential toroidal field with h = 0.1. (iv) Rigid toroidal field with k = 1. (v) Rigid toroidal field with k = 2. 


Table 2. Numerical results for each model. 


<? 

j2 

-'0 

h 

0 

MjT 

Smax (G) 

Smax / Save 

VAmax (cm/s) 

VC 

Rigid toroidal 

0.8 

7.12E-3 

- 

4.02E-2 

l.OOE-2 

3.67E16 

1.88 

4.14E8 

1.12E-5 

Type 2 differential toroidal 

0.8 

l.OOE-3 

0.5 

4.02E-2 

9.44E-3 

3.31E16 

1.75 

3.32E8 

1.12E-5 

0.8 

5.00E-4 

0.4 

4.02E-2 

8.25E-3 

3.09E16 

1.75 

2.94E8 

1.12E-5 

0.8 

2.00E-4 

0.3 

4.01E-2 

6.05E-3 

2.71E16 

1.79 

2.48E8 

1.12E-5 

0.8 

5.00E-5 

0.2 

4.01E-2 

2.91E-3 

2.08E16 

1.98 

1.86E8 

1.12E-5 

0.8 

3.50E-6 

0.1 

4.01E-2 

4.21E-4 

1.10E16 

2.76 

9.84E7 

1.12E-5 

0.8 

1.50E-7 

0.05 

4.00E-2 

2.83E-5 

4.56E15 

4.41 

4.07E7 

l.llE-5 
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Figure 2. Contours of p and the toroidal magnetic field (B^). The difference between the two adjacent contours is 1/10 of the maximum value of p or B^. (a) 
The contour of a p with non-magnetic field, (b) The contour of a Bip of rigid toroidal field with k = 1. (c) The contour of a Btp of type 1 differential toroidal 
field with h = 0.1. (d) The contour of a B^p of type 2 differential toroidal field with h = 0.1. 


■uj ~ 0.05 (h = 0.05) as seen in the left-hand panel of fig. [3] The 
profiles of Alfven velocities have the same tendency in the right- 
hand panel. The values of the ratio Bmax/Bave increases as the 
degree of the differential increases. The ratios are 1.75{h = 0.5), 
1, 98{h = 0.2), 2.76(/i = 0.1) and (4.41 h = 0.05). Thus, the dif¬ 
ferential toroidal magnetic field is localized and concentrated near 
the centre when the degree of the differential is large. 

However, the magnetic field energy ratios {M /T) of the dif¬ 
ferential toroidal field solutions are much smaller than that of the 
rigid one (Table [3. As fig. shows, the strength of the magnetic 
field and the Alfven velocity in each differential model near the lo¬ 
cal maximum point are almost the same as those of the rigid model 
at the same radius. However, the magnetic energies of the differen¬ 
tial models are smaller than that of the rigid model because of their 
localization and concentration near the point. Although the mag¬ 


netic energy ratio (M/T) of the rigid solution is ~ 1.0 x 10 
the energy ratios of type 2 differential solutions with small h are 
~ 2.91 X 10“® (h = 0.2), ~ 4.21 x 10"'' (h = 0.1), and 
~ 2.85 X 10“® (h = 0.05). The magnetic energy of the type 2 
differential toroidal solution with h = 0.05 is much smaller than 
that of the rigid toroidal model, even if the strength of the toroidal 
magnetic field is comparable to that of the rigid field at the local 
maximum point {■& ~ 0.05). The localisation and concentration of 
the toroidal magnetic field near the rotational axis can be realized 
by using type 2 differential toroidal fields. 
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Figure 3. Left: the profiles of B,p on the equatorial plane. Each curve denotes a rigid field with k = 1 (thin solid), and a type 2 differential field with /i = 0.5 
(thick dashed), h = 0.2 (thick dotted), /i = 0.1 (thick dot-dashed) and h = 0.05 (thick solid). Right: the profiles of Va on the equatorial plane. Each curve 
denotes a rigid field with k = 1 (thin solid), and a type 2 differential field with h = 0.5 (thick dashed), h = 0.2 (thick dotted), h = 0.1 (thick dot-dashed) 
and h = 0.05 (thick solid). 


4 DISCUSSION AND CONCLUDING REMARK S 

4.1 Suppression of the low-T /\W\ instability by low MjT 
magnetic fields 

The low -T/1 ml instabilit y is related to the corotation point inside 
the star ( lAnderssonll2003h . Since the wave pattern speed of an un- 
sTable mode is equal to the rotational velocity of the background 
star, the corotation resonance would be driven at the point. Recent 
time evolutions of linear analysis have shown that the f mode be- 
comes unsTable when the corot ation point appears inside the star 
dPassamonti & AnderssonI I 2 OI 5 I) . The low-T / |m| instability de - 
pends on the degree of the differential rotation dShibata et al.l2003h . 
but the location of the corotation point is not obvious. Linear anal¬ 
ysis has also illustrated that the location of the corotation point de¬ 
pends on_both_the_degree_of_thediffcrential rotation and the value 
of j3 dPassamonti & Anderssonll2015h . As the value of /3 decreases 
and the degree of the differential rotation increases, the location 
approaches the rotational axis. If the value of fi is 0.04 and the de¬ 
gree of the differential rotation is d = 0.3 as we have calculated in 
this paper, the corotation point will appear at to ~ 0.27 (see fig. 
6 in IPassamonti & Anderssonll2015l) . Therefore, the location of the 
corotation point will appear near the rotation axis in a differential 
rotating star. 

The suppression of the toroidal magnetic field may depend 
on the local strength of the toroidal magnetic field and Alfven ve¬ 
locity, since the slow-magnetosonic points are determined locally 
by the values of density and toroidal magnetic fields at the point 
dFu&Lal2011bh . As seen in figures [T] and the maximum point 
of the rigid toroidal magnetic field is located at to ~ 0.5 - 0.6. 
These values are larger than the location of the corotation point 
(to ~ 0.27). On the other hand, the maximum point of the type 
2 differential toroidal magnetic field is located af fhe inner region 
of the star. The maximum point is located at to ~ 0.2 when the 
degree of the differential is h = 0.2. The location approaches the 
rotation axis as the degree of the differential toroidal magnetic field 
increases (fig.[3]l. Since the differential toroidal magnetic field rep¬ 
resents the magnetic winding by the differential rotation very well 
(fig-[D, the profile of the toroidal magnetic field would become that 
of the type 2 differential toroidal magnetic fields if the toroidal field 
is wound up from the initial poloidal magnetic field by the differen¬ 


tial rotation. If the star has a large differential rotation, the degree 
of the differential toroidal magnetic field would also become large 
because the toroidal field is developed by the magnetic winding. As 
a result, both the maximum point of the toroidal magnetic field and 
the corotation point approach the rotational axis. We assume that 
both points appear near the centre when the star has a large differ¬ 
ential rotation. Such a differential toroidal field would suppress the 
low T /\W\ instability easily, even if the total magnetic field energy 
ratio (A4 jT) were much smaller than that of the rigid toroidal field, 
because the strength of the local toroidal magnetic field and Alfven 
velocity are comparable to those of the rigid toroidal field (fig. [3}. 

Such low-Af/T suppression mighf be realized by th e type 
2 differential toroidal field. Since iMuhlberger et alj ( 1201 4l) calcu¬ 
lated both rotation and magnetic winding simultaneously, the dis¬ 
tribution of their toroidal magnetic field would be type 2 differen¬ 
tial rather than rigid. The suppression of the low-T /\W\ instability 
might occur in a differential magnetic field configuration even if the 
magnet ic energy ratio is very small. On the other hand, IFu & Lail 
( I 2 OI Ibl) adopted a simple rigid toroidal magnetic field solution as 
a background model for their linear analysis. The suppression of 
the low-T /\W\ instability might need larger toroidal magnetic field 
energy than is available from a differential field because the maxi¬ 
mum point of the toroidal field is located at the outer region of the 
corotation point. 

The suppression would depend on the degree of the differen¬ 
tial toroidal field as the low-T /\W\ instability depends on the de¬ 
gree of the differential rotation. In order to study the low-Ad/T 
suppression systematically, in the future, we will perform linear 
analysis using the differential toroidal field obfained here. Such 
work may explain the difference in the energy ratio (M/T) be- 
tween a linear analysis with a rigid toroidal field dFu&LafcOllbl) 
and a numerical simu lation with a differentia l toroidal field made 
by magnetic winding dMuhlberger et al.l2014h . 

4.2 Concluding remarks 

We investigated new magnetized equilibria with differential rota¬ 
tion and differential toroidal fields in this paper. We developed 
new functional forms for differential toroidal fields in analogy to 
the rotation law of the differential rotation and obtained equilib- 
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rium states self-consistently by using a self-consistent numerical 
scheme. 

The dijferential toroidal magnetic fields have two functional 
forms. The functional form of the type 1 dijferential toroidal field 
is similar to the u-constant differential rotation law. The profile of 
the toroidal field is a gradual slope (fig.[I]l. Since the toroidal mag¬ 
netic field becomes constant when the star is incompressible, the 
type 1 dijferential toroidal magnetic field is analogous to a Bip- 
constant toroidal field law. The type 2 differential toroidal field is 
similar to the y-constant rotational law. Such differential toroidal 
fields represent the magnetic winding by the y-constant differen¬ 
tial rotation very well. The profiles of the toroidal magnetic fields 
represent the magnetic winding by the differential rotation (fig.|T]l. 
The type 2 dijferential toroidal field is localized and concentrated 
near the rotational axis when the degree of the differential toroidal 
magnetic field is large (fig.[2j. 

The low-r/|VK| instability would be suppressed by the con¬ 
centrated and localized dijferential toroidal fields more efficiently 
than a rigid toroidal field. The maximum strength of the differen¬ 
tial toroidal magnetic field can be locally comparable to that of a 
rigid toroidal magnetic field, even if the total magnetic energy of 
the differential field is much smaller than that of a rigid one (fig. 
O. Since the lo'x-T/\W\ instability and the suppression hy the 
toroidal magnetic field would depend on the strength of the local 
rotational velocity and the toroidal magnetic field, the localized and 
concentrated toroidal magnetic field is efficient for the suppression 
even if the total magnetic energy is small. Such low-AT /T suppres¬ 
sion would be realized by the differential toroidal field wound up 
from the initial poloidal magnetic field by differential rotation. The 
dijferential toroidal magnetic field plays a key role in explaining 
the difference in the energ y ratio (AT /T) be tween a linear analysis 
with a rigid toroidal field iFu & Laill20lT^ and a numerical sim¬ 
ulation with a differential toroidal field wound up from the initial 
poloi dal magnetic field by differential rotation jMuhlberger et al.l 
|2014 . 
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APPENDIX A: ACCURACY VERIFICATION 


To check the numerical accuracy, we computed a relative value of 
the Virial relation as follows: 


VC = 


\2T + W ^■m + M\ 
\W\ 


(Al) 


Our numerical domain is defined as 0 < f < 1 in the radial 
direction and 0 < 6 < n 12 m the angular direction. Fig. lAll 
displays the value of VC as a function of the number of grid 
points in the r-direction. To plot this graph, we computed the so¬ 
lutions with j-constant differential rotation and a type 2 dijferen¬ 
tial toroidal field {h = 0.3) changing the number of grid points 
in the r-coordinate, bu t fixing the number o f grid points in the 6- 
direction as ng = 513 jFuiisawa et al.|[20T^ . The value of VC de¬ 
creases in proportion to the s q uare inverse of the n umber of grid 
points jUander & Jones! l2009l : iFuiisawa et al.l l2012h . because we 
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log(N,) 


Figure Al. The Virial quantity VC, plotted against the number of grid 
points in the r-direction. 


adopted Simpson’s scheme as an integral method in our numeri¬ 
cal code. The typical value of VC is smaller than 10“^ — 10“® 
when the mesh numbers were Nr = 513 and Ne = 513. This 
val ue is sufficiently small and the system converges very well 
(cf. Hachisul 19S6^ Tomimura & Eriguchill200^ : iLtmder & JonesI 
1 20091 : Fujisawa & Eriguchill2013h . Thus, we fixed Nr = 513 and 
No = 513 in all calculations in this paper. 
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